function [bgRegionStd,bgRegionMean]=ODimageMeanStd(odimage, bgROI)
    %cal standard devation of background
%     bgROI=[667 160 300 300];
    bgRegion=odimage(bgROI(2):bgROI(2)+bgROI(4),bgROI(1):bgROI(1)+bgROI(3));
%     odimage(bgROI(2):bgROI(2)+bgROI(4),bgROI(1):bgROI(1)+bgROI(3))=10;
%     figure(10);imagesc(odimage);
    bgRegionStd=std(std(bgRegion));
    bgRegionMean=mean(mean(bgRegion));
%     rectangle('position',[bgROI(1) bgROI(2) bgROI(3) bgROI(4)],'Curvature',[1,0],  'EdgeColor',[1 0 0]);